Entanglement and the Interplay between Staggered Fields and Couplings 
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We investigate how the interplay between a staggered magnetic field and staggered coupling 
strength affects both ground state and thermal entanglement. Upon analytically calculating ther- 
modynamic quantities and the correlation functions for such a system, we consider both the global 
Meyer- Wallach measure of entanglement and the concurrence between pairs of spins. We discover 
two quantum phase transitions present in the model and show that the quantum phase transitions 
are reflected in the behaviour of the entanglement at zero temperature. We discover that increasing 
the alternating field and alternating coupling strength can actually increase the amount of entan- 
glement present at both zero temperature and for thermal states of the system. 
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I. INTRODUCTION 

Entanglement is an intriguing concept in quantum in- 
formation, and a resource used in many quantum com- 
putation schemes P. The Heisenberg coupling in spin 
chains has been shown to allow universal quantum com- 
putation Spin chains are also good candidates for 
quantum wires (3|. Thus it is an important task to ex- 
plore how the amount of entanglement in such systems 
changes under different conditions, to discover whether 
it can be enhanced and what causes its destruction. 

Many-body entanglement [l[ is difficult to quantify, so 
entanglement measures are often restricted to the pure, 
zero temperature case, for example von Neumann en- 
tropy Q or the Meyer- Wallach measure @ , or to a small 
number of possibly mixed qubits such as concurrence Q . 
On the other hand, it is, in general, hard to quantitatively 
study the entanglement of thermal states due to the ab- 
sence of efficient separability criterion for mixed states 
in many-body systems. An alternative approach is to 
use an entanglement witness which detects rather than 
measures entanglement. In particular, a thermodynamic 
entanglement witness uses thermodynamic quantities de- 
rived from the partition function of the system to detect 
entanglement 0, [1]. 

In spin chains, ground state entanglement is often in- 
vestigated in the context of quantum phase transitions 
(QPTs). A QPT is a sudden change in the properties of 
the ground state when a parameter of the Hamiltonian 
such as a magnetic field is varied. Since it is expected 
that by investigating QPTs, dramatic changes to physi- 
cal quantities at very low temperature would be revealed, 
QPTs have been intensely studied in spin systems Q . In 
Ref. [loj . it was shown that in general, a singularity oc- 
curs in the ground state entanglement at quantum critical 
points (QCPs). 

In this paper, we investigate the entanglement of a spin 
system in a non-uniform magnetic field with non-uniform 
coupling constants. Certain solid state systems such as 
Copper Bcnzoate have a non-uniform magnetic field , 



caused by an inhomogeneous Zeeman coupling. Similarly, 
examples exist for a non-uniform coupling strength |12l . 
[l3| . Some properties of these materials can be captured 
using the staggered spin chain we study in this paper. 

Entanglement in a staggered magnetic field has been 
studied previously using single-site entropy and an en- 
tanglement witness [i^ . The effect of such a staggered 
field on the dynamics and on the high-fidelity transfer 
of entanglement has also been considered [15|; it was 
found that the staggered field is almost as efficient as 
the uniform case. On the other hand, entanglement in 
a spin chain with staggered coupling and magnetic field 
has been considered only numerically for chains of finite 
length [11. 

Here, we show how to calculate two measures of en- 
tanglement analytically for such spin chains with infi- 
nite length. We first calculate thermodynamic proper- 
ties and finite temperature correlation functions of the 
spin chain. Using these results, we show that the sys- 
tem exhibits QPTs at zero temperature induced by the 
staggered fields. We investigate global entanglement of 
the ground state using the Meyer- Wallach measure, a 
measure of multipartite pure states based on bipartite 
entanglement. Wc then investigate both zero and finite 
temperature entanglement between two spins using the 
concurrence. We find that global entanglement in gen- 
eral increases with an alternating coupling constant. At 
zero temperature, the concurrence first increases, reaches 
a maximum, then decreases with an alternating coupling 
constant. We also find that for certain values of magnetic 
field, both the Meyer- Wallach measure and the zero tem- 
perature concurrence can be increased by the alternating 
magnetic field. At finite temperature, the concurrence 
can again be increased by certain values of the staggered 
fields and couplings. Further, we examine the entangle- 
ment at finite temperature using an entanglement witness 
in view of searching for thermal multipartite entangle- 
ment. 

This paper is organized as follows. In Section |lll 
wc present the Hamiltonian and give its diagonal form. 
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which allows us to calculate thermodynamical quantities 
of the system. We compute the correlation functions in 
Section lnil In SectionlTvl wc investigate the ground state 
and show that the Hamiltonian exhibits QPTs. Then, we 
investigate entanglement at zero temperature and at fi- 
nite temperature in Section |Vl We present concluding 
remarks in Section [VII 



II. THE HAMILTONIAN AND ITS 
THERMODYNAMIC PROPERTIES 

We consider the thermodynamic limit of the staggered 
Hamiltonian 



The partition fimction, Z = tr(e ^^), of this system 
can be written Z = Oa^i and thus In Z ~ J2k=i 



where Zk 
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is found from the fcth term in the Hamiltonian sum. Ex- 
plicitly taking the thermodynamic limit, wc find 

\nZ = — / dq\n [4 cosh (/3A+) cosh (^A")] (5) 
Jo 

where = B ± \/ cos^ q + b'^ + sin^ q. From this, 
other thermodynamic quantities such as the internal en- 
ergy, [/ = — ^ In Z, can be calculated: 
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(1) 



where Ji = J + e^^'-j is the staggered coupling strength 
and Bi = B + e^'^^h is the staggered magnetic field. We 
refer to j as the alternating coupling strength and b as 
the alternating magnetic field. This spin chain can be 
diagonalised [131 using a Jordan- Wigner transformation, 
'^i = rim^i ® + '''■f )/2, with anti-commutation 
relations {a;, = and {a[, a„} = followed by a 

Fourier transform, ai = N-^/'^Y^k'^ke^^''''^'^ ■ ^^ext the 
Hamiltonian must be rewritten as a sum to N /2: 



N 



r ^ [A+ tanh (/3A+) + A" tanh (/JA")] 
Jo ^'^ 



The magnetisation, M = InZ, is 



(6) 



^ [tanh (/3A+)+ tanh (/3A-)], (7) 
and the staggered magnetisation , Mg = \nZ, is 



Ms r dq b [tanh (/3A+) - tanh {l3A-)] 



N Jo 27r ^ J2 cos^ q + b'^ +f sin^ q 



(8) 



N/2 

^ = J2{t^k4'ik+ f4dl+N/2dk+N/2 (2) 
k=l 

+ v+dldk+N/2 + v^dl^^f^Jk-'2^Bl^ 

where /i^ = [2B ± 2Jcos(27rA;/Af)] and = [2b ± 
2«j sin(27rfc/A^)]. This can now be diagonalised using a 
canonical transformation, 

dk = Oik cos 9k + Pk sin Ok (3) 
dk+N/2 = -akiiinek+ Ihco^dk 

where Ok is determined by — Jcos(27r/c/A'') sin20fc + 
bcos26k = ij sin(27rfc/Af). Thus the diagonal form of 
the Hamiltonian is 

N/2 

^ = E {^t^^k + XkPlf^k - 251) (4) 

k=l 

where Xf = 2B ± 2^J^ cos^ (2ife) + 52 + ^-2 31^2 ^2^^^ 

The anti-commutation relations are now {aj,a„i} = 

{PjiPm} = (5;,™ and {a;,Q;„,} = {/3/,/3,„} = {a|',^„} = 
{a/, /3m} = 0. In all figures in the paper, we will fix J to 
1. 



III. CORRELATION FUNCTIONS 

We will use the correlation functions of the sys- 
tem to calculate each measure of entanglement. Since 
[H, erf] = 0, the only non-zero correlation functions of 
this spin chain are {crfcri+fj+<^fcri+ft), {^i<^i+r) and (erf). 
The Hamiltonian is scmi-translationally invariant in that 
all odd and all even sites can be considered identical in 
the thermodynamic limit. Due to this semi-translational 
invariance, when R is even, (erf) — (crf^^). 

We calculate each of the correlation functions following 
the method in Wc find these for any R. Generally, 
{^i^i+R + = 2(aJ 11^=1+1(1 - 2al^a,n)ai+R + 

«l+Hnm=l+i(l-2aJ„a,„)ai), (CTfcrf_^H> = {{'^-'^a\ai){l- 
2al_^_j^ai+R)) and (erf) = (1 - 2a[a;). 

The auto-correlation functions of this model have been 
found analytically in the infinite temperature limit [l9| . 
and the time dependent (erf (i)CTf_^^) has been found, also 
analytically, for arbitrary temperature [20| . However, the 
general equilibrium correlation functions have not been 
calculated previously. 

Since the spin chain is a free fermion model and scmi- 
translationally invariant. Wick's theorem can be applied 
to the total correlation functions. Thus we can rewrite 
each of the above equations in terms of two point corre- 
lation functions. For example, the zz correlation func- 
tion is ((1 - 2alai){l - 2al_^^^ai+B)) = {<j!){crf^j^) - Gl„ 
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where we have defined Gj.ij 
We use the notation Gi^n = 
(cr^)° + e"'(cr^)'' where G?; 



Gfj. and (ct, 



and G! 



'■'ajai+R - aiaj^fi) (similarly for 



(af )). Note that we have treated (af) separately to Gi^n 
despite the fact that R = should give us (erf). The 
reason for this will become clear below. 

Using the Jordan- Wigncr transformation and Fourier 
transform, and then summing to N/2 (the additional 
canonical transformation is unnecessary), we find we can 
write 



N/2 



k=l 
N/2 

E[(4 



1 - 2dldk) 



(1 - 2dl , „ d 



2 

'n 



k=l 



(9) 



for (erf) and 



B=0.5,j=0.3 



B=l.S,j=0.5 




rr 



(b) 



FIG. 1: Second derivative of the ground energy with respect 
to b. The figures (a) and (b) are for {B,j) = (0.5,0.3) 
and {B,j) — (1.5,0.5), respectively. It is observed that the 
second derivative diverges at the QCPs B — \/ + b'^ and 

B = ,/pT¥. 



We note that for R even, G|j. = 0. For odd R, we in- 
stead compare the correlation function form of Gi.r to 
Gi^i, again noticing they are similar. Using the thermo- 
dynamic form of G/_i, we determine the thermodynamic 
form of G; jj: 



N/2 



1^ R 



-y 

k=l 



(2'KkR\ 



\ N ) 



[(1 



(10) 



r dq , Jcosg[tanh(/3A+) -tanh(/3A-)l 
Gr = - — cos(gi?)- 



Gf 



2tt 



62 



•9 • 2 

J ^ sm q 



COS"" q 

^ dq j sin q [tanh(/3A~'" 

27r yJ^ cos^ q + b'^ + p s\t? q 



tanh(/3A-)] 



Gf 



2^^'^ 



k=l 



2TrkR 
N 



(44-i-f ) 



A-kR. 



{dl^N^k) 



The differences between odd and even R are due to the 
presence of the e*'^^ term in the correlation function 
forms of G;.ij. 



for Gi,R. We know the thermodynamic forms of (trf )° = 
m (Eq. El), (erf)" = m, (Eq. [5]), and also of G/,i = 
— {afafj^^ + a^a^^-^)/2 which can be calculated by dif- 
ferentiating the partition function with respect to the 
couphng strengths, G°j = — -^^InZ and Gf^ = 

Thus we have 



Gl 
G\ 



r dq Jcos2 q [tanh(/3A+) - tanh(/3A-)] 

/o 27r ^ J2 (,Qg2 q _|_ ^2 _j_ j2 gjjj2 ^ 

r dq j sin^ q [tanh(/3A+) - tanh(/3A-)] 
Jo 27r ^ J2 cos2 q + 62 + j2 sin2 ^ 



The thermodynamic expressions for Gi_r are found di- 
rectly from the above equations. For even i?, we com- 
pare the correlation function form of Gi^r to that of (erf), 
noticing that they are similar. Since we know the ther- 
modynamic form of (erf), we can also determine the ther- 
modynamic form of Gi,r: 



Ltd 



Gs 
R — 



dq 
2^ 



.s(gi?) [tanh(/3A+) +tanh(/3A-)] (11) 



dq bsin{qR) [tanh(/3A+) - tanh(/3A-)] 



2tt 



\/ J2 cos2 g + 62 + j2 gjj^2 ^ 



IV. PROPERTIES OF THE GROUND STATE 

In this section, we investigate the properties of the 
ground state. Without the alternating coupling strength 
and the alternating magnetic field, the Hamiltonian is 
referred to as an XX model with a transverse magnetic 
field. It is well known that the XX model has a second 
order QPT at i? = J [2l|. We investigate the ground 
state energy as a function of the alternating coupling 
strength and the alternating magnetic field, and see the 
QPTs induced by them. 

We first define Q as 

g-{Qe [0,7r] I A- <0}. (12) 
For simplicity, we also introduce two functions. 



e(g) 



J2 COs2 g + 62 



•9 • 2 



'52-62 



= arccos < 



J2 



(13) 
(14) 



By using the internal energy per site u given by Eq. ([SJ, 
the ground state energy per site is obtained by taking the 
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(a) B=0.6 (b) B=1 




FIG. 2: Magnetization at zero temperature. The figures (a), 
(b) and (c) sliow B — 0.6, B — 1 and B = 1.4, respec- 
tively. The magnetization changes non-smoot hly at th e quan- 
tum critical points B — y' J'^ + b'^ and B = + b^. 



limit such as 

eg = - lim / p- [A+ tanh(/?A+) + A" tanh(/3A-)] 

I3~^00 Jq ^ H H H 

A^Q 27r J^gQ 27r 

= -B + - [ [B-e{q)]dq. (15) 

In Appendix |X1 the region Q is analytically obtained and 
exact expressions of the ground state energy are shown. 
In order to visualize the QPTs, the second derivative of 
the ground state in terms of b for {B, j) = (0.5, 0.3) and 
{B,j) ~ (1.5,0.5) is given in Fig. [T] It is observed that 
the derivative diverges at the points B = yj^+l? and 
B = Vj+^i both of which lead to second order QPTs. 

The QPTs are observed more clearly by looking 
at the magnetic susceptibility at zero temperature. 
In Appendix \^ the magnetization at zero tempera- 
ture is explicitly calculated, which is a function of S. 
Since the derivative of S by i? and b includes a fac- 
tor 1/a/B2 - b'^ - fV-B'^ + J2 the derivative di- 
verges at the points B ~ \/ + b'^ and B — y^p + 6^, 
which implies QPTs. To demonstrate this, we also plot 
the magnetization at zero temperature in Fig. [21 We can 
clearly see that the magnetization of the ground state 
changes non-smoothly at both QCPs. 

We note that, when j ~ b = corresponding to the XX 
model with a transverse magnetic field, B = \\P+1}^ 
reduces to the QCP of the XX model, B = J. On the 
other hand, QPTs a,t B = \/j^+l>^ do not appear in the 
XX model. Hence this QCP can be regarded as being 
induced by the staggered nature of the spin chain. 



V. ENTANGLEMENT 

We study the entanglement properties of the spin chain 
in three different ways. First we define the Meyer- 
Wallach measure and the concurrence and discuss how 
to calculate each for the staggered Hamiltonian. Next we 
consider how both entanglement measures behave at zero 
temperature, and then how a finite temperature affects 
the concurrence. Finally we consider a thermodynamic 
entanglement witness in an attempt to detect thermal en- 
tanglement which the measures miss such as multipartite 
entanglement. 

A. Meyer- Wallach measure and Concurrence 

First, we show how to calculate the Meyer- Wallach 
measure and the concurrence from the thermodynamic 
quantities calculated in Section |TT] and correlation func- 
tions in Section Hill 

For a pure state |$) of an A^-spin system, the Meyer- 
Wallach measure is defined by 

2 ^ 

i?Mw(|$))=2--^Trp2, 

i=l 

where pi := Tr-^|$}($| is the reduced density matrix at 
the i-th spin 0, [23]. (The partial trace Tr^^ is taken 
over all degrees of freedom except the z-th spin.) The 
Meyer- Wallach measure takes values between and 1. 
The minimum is achieved if and only if the state is sep- 
arable while the maximum is given by states which are 
local unitary equivalent to the GHZ state. 

We calculate the entanglement of the ground state 
found by the Meyer- Wallach measure. Since the Hamil- 
tonian preserves the magnetization, i.e., [H,J2i'^i] — 
the reduced density matrix of a spin at site I, pi, has only 
diagonal elements such that pi = diag{i^t^|^, -!— ^|^-^} 
where (cf )g := {g\<7i\g) is the expectation value of af by 
the ground state \g). By substituting this and using the 
fact that the Hamiltonian is semi-translationally invari- 
ant, the Meyer- Wallach measure of the ground state is 
obtained as 

EMw{\g)) = 1 - lii'yLenfg + {<dd)l) ■ (16) 

where (cr|^en>9 = {g\(^li\g) and (cr^dd)s = (.9k2;+il5> for 
any I. The (cr^yen) and {cTodd) are obtained from the mag- 
netization m given by Eq. ([71) and given by Eq. ([H) 
such as 

{(^lven)g = lim [m + m,], (17) 
Kdd>s = lim ["^- "Is]. (18) 

The exact expressions of the ground state energy, the 
magnetization at zero temperature and the Meyer- 
Wallach measure are given in Appendix \^ Note that 
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the Meyer- Wallach measure is a measure of entanglement 
only for pure states, and thus is meaningful for investi- 
gation of entanglement in the ground state but not that 
of thermal states. 

The concurrence, C, between two spins @ is an entan- 
glement measure for both pure and mixed states, and can 
therefore be used at finite temperatures. It is given by 

C(p) = max{0, Ai - A2 - A3 - A4}, (19) 

where the A^s are the square roots of the eigenvalues of 
the matrix pp with p = (cr^ ® (jy)p*{ay (8) cr^), and sat- 
isfy Ai > A2 > A3 > A4. Again using [H^Ylii^i] = 0, 
the concurrence between two spins at sites I and / -I- i? is 
C{pi^i+r) = 2max{0, \z\ - ^/vy} where z = jiafcrf+n + 
afa,%),and^;2/=3^[(l + (afa,^«))2-((af) + (af+^))2]. 
Due to the semi-translational invariance of the Hamil- 
tonian, the concurrence is the same for all odd sites 
and for all even sites. The concurrence is zero if and 
only if the state of the two spins is separable, and is 
one when they are maximally entangled. Although we 
could calculate the concurrence for any R, we concen- 
trate on the nearest neighbour, Ci with R — 1, and 
the next nearest neighbour, C2 with i? = 2, concurrence 
since for large R the concurrence is infinitesimal. Us- 
ing {afaf^j{) = {(Tf){af^ji) - Gf_^, the nearest neighbour 
concurrence is 




(20) 

and the next nearest neighbour concurrence is 



C2 = max|0,|G,,iG,+i,i-G;,2('Tf+i)| (21) 

- ^\/(l + (-,X+2))'-(2(^f))'}, 

remembering that Gr is different for odd and for even 
values of R. 

When J = j, the total coupling strength between near- 
est neighbour sites for odd I is zero, while for even /, it 
is 2J. Thus, at any temperature, there is no entangle- 
ment between nearest neighbours for odd sites, and at 
zero temperature (and magnetic fields), the chain con- 
sists of A^/2 maximally entangled singlet states. This is 
an example of dimerisation. As a consequence of this, 
there is also no concurrence at J ~ j for both odd and 
even sites for any R > 1. 

We note that unlike the Meyer- Wallach measure, the 
concurrence is not directly related to the total amount 
of entanglement contained in a pure state. Hence, there 
is no guarantee that the Meyer- Wallach measure and the 
concurrence will behave similarly. For instance, for the 
GHZ state, the Meyer- Wallach measure is one but the 
concurrence between any two spins is zero. 



(a) j=o (b) j=0.5 




FIG. 3: The Meyer-Wallach measure as a function of B and 
b. 



B. Zero temperature 

In this section, we study entanglement at zero temper- 
ature using the Meyer-Wallach measure and the concur- 
rence as defined above. Figs. |3] and [S] show the Meyer- 
Wallach measure. Figs. |4] and [6] the nearest neighbour 
(NN) concurrence for both odd and even sites, and Fig. [7] 
the next nearest neighbour (NNN) concurrence. 



1. Quantum phase transitions 

We first discuss the Meyer-Wallach measure of the 
ground state. See Appendix |A] for the detailed calcu- 
lation. In Figs. [3] and O it is observed that the Meyer- 
Wallach measure changes non-smoothly at the QCPs. 

Conversely, considering the plots for concurrence at 
zero temperature, the quantum phase transitions present 
in the staggered model are not always evident. In par- 
ticular, when plotting B against b for odd sites in Fig. |4j 
only the curve B ~ can be seen in (b), while 

only B = \J ]^ + \P- can be seen in (c). However, changes 
in the concurrence can be observed at both QCPs for 
even I in the same figure, and for both odd and even 
sites in Fig. [5] where we plot B against 

These results are consistent with the results in 
Ref. [l^l, where it is shown that entanglement of the 
ground state behaves singularly around QCPs in general. 



2. Effects of j and b on entanglement 

Next, we discuss the effects of b and j on the Meyer- 
Wallach measure and the concurrence. 

First, we consider the effect of b on the entanglement, 
shown in Figs. [3] and 21 How the concurrence behaves 
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(b) j=0.5 I odd 



(a) b=0 odd I 



(b) b=0.5 odd I 
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FIG. 4: The nearest neighbour concurrence as a function of 
B and 6 as T — >■ 0. When j = 0, (a), the concurrence for odd 
and for even sites is identical. 



(a) b=0 



(b) b=0.5 
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(f) b=0.5 even I 
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FIG. 5: The Meyer-Wallach measure as a function of B and 



2 3 4 5 
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FIG. 6: The nearest neighbour concurrence as a function of 
B and j as T — >■ 0. (i) plots the concurrence against j for 
I even, B — 0.5 to show that the concurrence decreases for 
high enough j for 6 = (solid red line), b — 0.5 (dashed green 
line), 6=1 (dotted purple line) and 6 = 1.5 (dot-dashed blue 
line) . 



in the presence of the alternating fields is highly de- 
pendent on whether the site is odd or even. In gen- 
eral, due to the larger coupling from an even site to 
an odd site, both nearest (NN) and next nearest neigh- 
bour (NNN) concurrence is higher, with the opposite be- 
ing true from an odd to an even site. Entanglement 
for the Meyer-Wallach measure and odd site NN con- 
currence remains large only when B is between the two 



QPTs, i.e. y^p + 62 < B < y/,P + b'^ when j < J or 
y/ + b'^ < B < \J p + &2 when J < j. For both mea- 
sures, the maximum entanglement is at B ^ 6 for large 
magnetic fields. These results are understandable from 
the fact that a large magnetic field aligns spins in the 
same direction leading to a separable state. When 6 and 
B are large, entanglement can be large only for B ^ b 
since the magnetic field on odd sites is canceled in such 
cases. 
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(a) j=0 odd I 



(b) j=0.5 odd I 
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of maximally entangled state. 

Another interesting feature common to both the 
Meyer- Wallaeh measure and the concurrence, demon- 
strated in Figs.[3]-|ni is that below the region of the QCPs, 
the entanglement is constant as B varies. That is, the 
magnetic field B is not a dominant parameter for entan- 
glement below the QCPs. The QPTs are often intuitively 
understood as occurring due to the balance between the 
strength of the coupling constants and that of the mag- 
netic fields. As such, the dominant parameters of this 
system are the coupling constants (the magnetic fields) 
below (above) the QCPs in general. Our results sup- 
port this intuition from the viewpoint of entanglement 
in the sense that the magnetic field B does not change 
the entanglement below the QCPs. On the other hand, 
the entanglement is sensitive to the change of the al- 
ternating magnetic field b even below the QCPs, which 
demonstrates the difference between B and b. 

Finally, the NNN concurrence is reduced compared 
to NN concurrence as expected, but remains reasonably 
high, especially for even sites where a non-zero value of b 
increases the entanglement. Further, increasing 6 allows 
a spin chain with larger values of B to be entangled. 



FIG. 7: Next nearest neighbour concurrence as T — 0. 
currence is zero when j = 1 as discussed in the text. 



Con- 



C. Finite temperature 



The even site NN concurrence does not follow this pat- 
tern, and instead a larger amount of entanglement tends 
to be present when B < \J ■\- 6^ for j < J or when 
B < V.P + b^ for J < j. 

On the other hand, the Meyer- Wallaeh measure and 
the concurrence vary differently with the alternating cou- 
pling constant j as shown in Figs. [S] and [5] where the 
entanglement is plotted as a function of j. The Meyer- 
Wallach measure is an increasing function with j except 
in the vicinity of the QCPs while the concurrence is in 
general a decreasing function of j. Fig. [5] (i) shows that 
for each b, there is a non-zero value of j for which the 
concurrence is the maximum possible. Since the Meyer- 
Wallach measure is a measure of the entanglement shared 
in the whole spin chain and the concurrence measures the 
entanglement between two spins, we can conclude from 
these results that as the alternating coupling constant j 
increases, the amount of entanglement shared amongst 
all spins increases. Such global entanglement is not lo- 
cally detected, in the sense that the entanglement of the 
reduced two spin state is small. 

The concurrence is closely related to the entangle- 
ment of formation. In Fig. [51 for odd sites, when 
B > + 5^, increasing j can increase concurrence, 
and for even sites, when B < \J p -\- increasing j in- 
creases the concurrence until a maximal value is reached, 
after which the concurrence decreases again. Thus a low 
but non-zero value of j can be beneficial to the extraction 



We next investigate the entanglement properties of the 
thermal states of the Hamiltonian using the concurrence. 

Fig.[8]plots the nearest neighbour concurrence for both 
odd and even sites for varying temperature and alternat- 
ing coupling strength, j. The figures show that increasing 
j allows the spin chain to be entangled at higher temper- 
atures, and that increasing both b and B can enlarge the 
region of entanglement. That is, the spin chain is entan- 
gled for more values of T and j for higher b and B. This 
is true for even as well as odd sites. Next nearest neigh- 
bour concurrence can be seen in Fig. [3] where we again 
see that increasing the magnetic fields can be beneficial 
to entanglement. As discussed previously, there is no 
entanglement at j = I for odd sites at any temperature. 

Increasing temperature has the effect of mixing energy 
levels, something which has the ability to either increase 
or decrease entanglement, though a high enough temper- 
ature will destroy entanglement. Increasing the alternat- 
ing coupling strength counteracts this to some extent, 
though a high enough temperature will still destroy the 
concurrence. 

We note that increasing b allows larger values of con- 
currence at higher temperatures at lower, more accessible 
values of j as demonstrated in Fig. |S] (b) and (f). Thus 
it is the combination of alternating fields b and j that 
allow the spin chain to be entangled at higher temper- 
atures. Increasing B has a similar effect, though to a 
lesser extent. 
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(a) b=0 B=0 odd I 



(b) b=5 B=0 odd I 
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(f) b=5 8=0 even I 




(h) b=5 B=5 even I 




FIG. 8: Nearest neighbour concurrence for odd and even sites. 



(a) b=0 B=5 



(b) b=5 8=5 even I 




FIG. 9: Next nearest neighbour concurrence. When 6 = 
and B — 5, (a), the plot for odd and even sites is identical. 
No next nearest neighbour concurrence is found for the other 
values of B and b shown for nearest neighbour concurrence. 



D. An entanglement witness 

In order to detect rather than measure entanglement 
in this system, we use an entanglement witness based on 
the expectation value of the Hamiltonian: 




FIG. 10: We witness the entanglement at different values of 
the alternating coupling strength, j for j — purple, j = 0.5 
blue, j = 1 red, j = 1.5 green, from inside to out. 



4\U + BM + bMs 
N{\J-j\ + \J + j\) 



< 1, 



(22) 



where U is the internal energy (Eq. |6]) M is the 
magnetisation (Eq. [7]) and Ms is the staggered mag- 
netisation (Eq. [S]). The bound is found similarly to 
the usual method 0, Q- Rearranging the expectation 
value of the Hamiltonian gives 2\U + BM + 6M.J = 



'I "i+i 



BM 
+ \J 



bMs\/N < |(J - j)||(afaf+i + afaf^,, 
j\\((jf<jf^i +<Tfaf, i)li.R„Ri,. Next, the bound for both the 



The absolute sign allows us to write 4|[/ 

(o-ro-r+l+crfcrf+i)|i,e.„en- " 

odd and even / for pure product states can be found using 
the Cauchy-Schwarz inequality, and the definition of the 
density matrix giving \{o'f(jf_^i + afaf_^^-^)\ < 1. Due to 
the convexity of the set of separable states, this bound is 
also true for all separable states while an entangled state 
can violate this bound. 

Fig. [TU] demonstrates that increasing the alternating 
coupling strength increases the region of entanglement 
detected by the witness. That is, at larger j, entangle- 
ment is detected for higher values oi B, b and T than is 
possible at smaller j. This trend persists even at very 
high values of j. In addition, the alternating magnetic 
field increases the maximum values of B for which entan- 
glement is detected. However, overall, increasing i?, b or 
T enough (except at B ^ 6 as discussed in Section rVB[) 
will destroy entanglement either by causing the spins to 
align with the magnetic field, or via the mixing of energy 
levels as the temperature is raised. 

The entanglement witness generally complements the 
results of the entanglement measures, and allows for the 
possibility of detecting multipartite entanglement that 
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cannot be measured by them. However, for our Hamilto- 
nian, comparing Fig. [TUlto Figs. H] and 111 it can be seen 
that this witness does not detect any extra entangled re- 
gions compared to the concurrence. 



1. Ground energy 

First, we show the exact expressions of the ground en- 
ergy tg given by Eq. (jisp . 



VI. CONCLUDING REMARKS 



We have found that the introduction of an alternating 
coupling strength and alternating magnetic field into the 
usual XX spin chain in a uniform magnetic field can, 
for certain values of the parameters, increase both the 
amount and region of entanglement quantified by either 
the Meyer- Wallach measure or the concurrence. This is 
the case for both zero and finite temperatures. We have 
demonstrated that two quantum phase transitions exist 
in this system, signs of which are evident in both entan- 
glement measures. In addition, we have calculated an 
entanglement witness which detects entanglement within 
a region which agrees with the measures of entanglement 
we consider. 

It would be interesting to calculate the finite temper- 
ature effects of the quantum phase transitions in this 
model. Determining the effect on entanglement of in- 
creasing the period of the staggered parameters would 
also be an interesting extension to this work. For exam- 
ple, by varying the magnetic field and coupling strength 
over three sites /, Z -I- 1 and / + 2 rather than the two 
considered here. However, this may not be possible to 
do analytically. 
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In this case, Q is given by 

0= (0,S)U(7r-S,7r). (A5) 
It is straightforward to calculate the ground energy; 

^ ^/(fS-l)B-f /o-e(g)d<z, forSeBi 



for B e B'^. 



; (A6) 



bmce A- ^ 2B - 2VJ2 + fe2^ the ground energy is 



obtained as 



-VJ2+F, for BeBi 



for B e B'.. 



(A7) 



In this case, Q is given by 

Q = (S,^-5). 
Then, the ground energy is calculated as 

\~!^'^ ^{l)d<l^ for Be Si 

eg = \-l^B-l /|/' Q{q)dq, for B & B^ 
y-B, for Be 63. 

2. Magnetization at zero-temperature 



(A8) 



(A9) 



Appendix A: Calculations of quantities at zero 
temperature 

Here, we give exact expressions of the ground energy, 
the magnetization at zero temperature and the Meyer- 
Wallach measure of the ground state. For simplicity, we 
define regions such as 



B, := [0, + 62), 

52 := W.P + h\ /fT62), 

53 := [/7^T6^,oo), 



63 := 



[v/JM^,oo) 



(Al) 
(A2) 
(A3) 
(A4) 



We show the magnetization at zero temperature per 
site nig. The magnetization nig is directly obtained from 
Eq. ^ such as 

TO = lim / ^[tanh(/3A+) +tanh(/3A")] (AlO) 
/3^ooJq 27r 



dq 



(All) 



By substituting Q, the magnetization is obtained as fol- 
lows; for j < J 



1-^S, for Be 61 
1, for B eB'. 



(A12) 



3' 
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for j = J, 



For j = J, 



0, for B eBi 

1, for B eB'. 



3' 



(A13) 



and, for j > J, 



0, for B e Si 
ruo = < f S, for B GB2 

1, for B e 63. 



(A14) 



MW 



J2 + b2 



for B G Si 
for B e 6^. 



(A16) 



3. Meyer- Wallach measure of the ground state 



Finally, for j > J, 



Here, we give the exact expressions of the Meyer- 
Wallach measure of the ground state, Emw-, given by 
Eq. dUl). For j < J, 



MW 




(v/o^efe^?)^ for Be Si 
for B e B'3. 



(A15) 



E, 



MW 



.0 for Be 63. 



(A17) 
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